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We perform numerical renormalization group (NRG) as well as analytical calculations for the two- 
channel Kondo model to obtain the dependence of the Kondo temperature Tk on the dimensionless 
(bare) spin exchange coupling g over the complete parameter range from g <C 1 to g ^> 1. We show 
that there exists a duality between the regimes of small and large coupling. It is unique for the 
two-channel model and enables a mapping between the strong and the weak coupling cases via the 
identification g <-> 3/(2$), implying an exponential dependence of Tk on 1/g and g, respectively, 
in the two regimes. This agrees quantitatively with our NRG calculations where we extract Tk (<?) 
over the complete parameter range and obtain a nonmonotonic TK(g) dependence, strongly peaked 
at the 2CK fixed point coupling g* . These results may be relevant for resolving the long-standing 
puzzle within the 2CK interpretation of certain random defect systems, why no broad distribution 
of Tk is observed in those systems. 

PACS numbers: 72.10.Fk, 72.15.Qm 



I. INTRODUCTION 

The Kondo effect^ is a paradigm for strong electronic 
correlations in metals, induced by resonant quantum spin 
scattering of electrons at the Fermi energy from local de- 
fects with spin 5*. The generalization of the problem 
to the case of M equivalent conduction electron chan- 
nels, the multi-channel Kondo problem, has attracted 
much attention ever since it has been introduced by 
Nozieres and Blandin 2 in 1980. While for a channel 
number M = 2S the impurity spin is exactly compen- 
sated by the conduction electron spins below the Kondo 
temperature Tk, corresponding to a spin singlet strong 
coupling fixed point with Fermi liquid behavior—, they 
showed that for M > 2S both the weak and the strong 
coupling fixed points are unstable, and hence a stable in- 
termediate coupling fixed point was conjectured. It cor- 
responds to an overcompensation of the impurity spin 
at low temperatures due to the simultaneous screening 
by each channel, implying a nonvanishing zero-point en- 
tropy and non-Fermi liquid behavior. In the following we 
will focus our discussion on the spin S = 1/2 two-channel 
Kondo (2CK) effect. The anomalous behavior of various 
thermodynamic quantities near the 2CK fixed point has 
been worked out theoretically using the Bethe ansatz^^, 
a Majorana fermion representation of the problem^ and 
conformal field theoryj 7 -^ Early on, the two-level-system 
(TLS) model of atomic defects embedded in a metallic 
host was put forward by Zawadowski and Vlada r 10 : 11 as 
a physical realization of 2CK defects, where the internal 
TLS degree of freedom takes the role of the Kondo spin 
(pseudospin) and the magnetic conduction electron spin 
serves as the conserved channel degree of freedom. How- 
ever, it was shown thereafter that, unfortunately, within 
this model the 2CK fixed point cannot be reached be- 
cause of the instability of the 2CK fixed point with re- 
spect to external perturbations: Within this model the 
TLS tunneling attempt frequency sets the band cutoff for 



the 2CK effect, since band electrons at higher energies in- 
stantaneously screen the tunneling defect without pseu- 
dospin flip. This turns out to prevent Tk to be greater 
than the tunneling-induced level splitting of the TLSJ^ 
It remains to be seen if this problem can be overcome 
by a recently proposed modified TLS model^ where the 
2CK fixed point may be stabilized by an additional res- 
onance enhancement of the conduction electron density 
of states (DOS). 

On the experimental side, signatures consistent with 
the 2CK effect have been observed in both, certain bulk 
heavy fermion compound a 14 ! 15 ' 16 and in mesoscopic de- 
fect structuresi 17 i 18 The existence of TLS fluctuators in 
nanoconstrictions has been established by various exper- 
iments ! 19 i 20 One of the best-studied case of 2CK signa- 
tures is perhaps the zero-bias conductance anomaly ob- 
served by Ralph et al. in nanoscopic point contacts of 
simple metalsi 17 i 18 ' 21 A scaling analysis of the differential 
conductance of these contact o 22 i 23 and systematic param- 
eter variations lend strong support to the 2CK hypoth- 
esis. However, the 2CK interpretation of these data has 
remained controversial 2 ^! 2 ^ due to the lack of an estab- 
lished microscopic model for the physical realization of 
the 2CK defects. See Ref. [26| for an alternative, statisti- 
cal explanation of the zero bias anomaly. Most recently, 
2CK behavior seems to have been realized by systemat- 
ically tuning a semiconductor quantum dot system into 
the 2CK regime^ 7 - as proposed theoretically in Ref. [28|. 

One of the problems with the 2CK interpretation of the 
anomalies in disordered, mesoscopic nanoconstrictions is 
the fact that within this interpretation these systems ex- 
hibit a sharp value of the Kondo temperature Tk , while 
one expects a broad distribution of the pseudospin flip 
coupling J due to the random nature of the 2CK defects. 
In fact, for single-channel Kondo impurities in nanocon- 
strictions the observed behavior— has consistently been 
explained^ in terms of a broad Tk distribution, induced 
by mesoscopic fluctuations of the local DOS. 
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In the present paper we make a contribution to the 
resolution of this puzzle. We compute the dependence of 
Tk{J) on J within the generic, symmetric 2CK model, 
covering the complete range from small to large J. Since 
the 2CK fixed point is at an intermediate coupling J*, 
one expects that for J = J* the 2CK regime extends in 



energy up to the band cutoff D 



11,31 



i.e. for the 2CK 



case Tk{J) should have a maximum at J — J* with 
Tk(J*) — D. In Section II we define the model and, 
following the ideas of Nozieres and Blandin, 2 establish 
a duality between the large J and the small J region 
which makes it possible to give analytical expressions for 
Tk(J) in both regimes. Details of this calculation can 
be seen in the Appendix. In addition, wc compute the 
complete dependence Tk{J) using the NRG, as explained 
in Section III. The results are presented in Section IV, 
which are in quantitative agreement with the analytic 
expressions of section II and indicate a strongly peaked 
dependence of Tk on J. The conclusions and possible 
consequences for the 2CK interpretation of anomalies in 
nanoconstrictions are drawn in Section V. 



II. DUALITY OF THE 2CK WEAK AND 
STRONG COUPLING REGIMES 

We consider the isotropic 2CK Hamiltonian, 



H 2C k = £ fc c L-< 



kaa 



where 



are the usual creation operators for electrons 



in channel number a = ±1 with momentum k and spin 



is the creation operator for an 



electron at the impurity site, a aa i the vector of Pauli ma- 
trices and S the impurity spin operator of size 1/2. The 
exchange coupling J > is taken to be antiferromagnetic. 
We define the dimcnsionless coupling g = poJ, where 
po = 1/2D is the DOS at the Fermi level. Throughout 
this paper, all energy scales and coupling constants are 
given in units of the band cutoff D. 

In the weak coupling regime, g <C 1, the crossover scale 
to the 2CK non-Fermi liquid behavior can be obtained 
by perturbative analysis in g. It is well known as the 
weak coupling Kondo temperature and reads, including 
subleading logarithmic corrections,- 3 " 



T^ C) ~ Ce -^+ln(2 S )+0( ff ) j g <<; ! 



(2) 



Turning now to the strong coupling regime, j> 1, it 
is convenient to consider the Hamiltonian in site repre- 
sentation, 



\ " 

kaa 



/ j ^^^kacr^kacr 



t ^ ] C iaa C jaa ' 
(i,j)acr 



where i is the site index and an infinite, one-dimensional 
lattice with a nearest neighbor hopping amplitude t is as- 
sumed without loss of generality. In the limit g —* oo the 



kinetic energy in the Hamiltonian Eq. ([T]) is negligible, 
and we have, 



„{sc) _ J _ V- f 3 
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(3) 



The mapping of the strong coupling regime of the 2CK 
model ([T]) onto a weak coupling problem proceedes in two 
steps. We first represent the Hamiltonian fl| in the basis 
of low-lying eigenstates of the strong coupling Hamilto- 
nian ((3]) , which will be of the type of a generalized Ander- 
son impurity model. Then we project this model in the 
low-energy regime onto an effective weak coupling 2CK 
model. 

The ground states of this strong coupling Hamiltonian 
© are 3-body states comprised of one electron in each 
of the two channels, located at the impurity site and an- 
tiferromagnetically coupled to the impurity spin. These 
states are easily calculated as 
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T/ __(2| 



T1U» 



Fhvac) (4) 



|^)=-L(2|m)-|m)-| W)) - f>ac}(5) 

and have the energy E = - J, H^\^> ni) ) = -J|*° a) ). 
In the Dirac ket notation above the thick arrow repre- 
sents the impurity spin, while the first and the third 
(thin) arrow describes the conduction electron spin in the 
a = +1 and a = — 1 channel, respectively. For later use 
we have also defined fermionic operators which create 
these states out of the vacuum \vac) (i.e. the free Fermi 
sea without impurity) . Note that the ground states can- 
not simply be product states of 2-particle singlets, but 
necessarily contain triplet admixtures, a frustration ef- 
fect implied by the quantum nature of the Hamiltoniam 
©. The degeneracy of the l^m) i s ^ ne reason wn Y 
the 2CK model remains nontrivial even in the strong 
coupling limit, in contrast to the single-channel Kondo 
model. The next excited eigenstates of Eq. © are the 
2-body singlet and triplet states \^ sma ), 
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V2 
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(I n o> - 1 ifr o» 
(I n o> + Ufto» 



(6) 

(7) 

(8) 
(9) 



and analogous definitions for the a = — 1 channel. In 
the above notation, s = 0, 1 denotes the total spin, 
m = 0, ±1 its z-component and a = ±1 the occupied 
conduction channel of the 2-body state. The energies of 
these states with respect to Eq. §5§ are E n0 a = — jJ and 
E\ ma — +j J, respectively. Switching on the hopping t 
removes an electron from the 3-body states Eqs. Q and 
puts it onto a site i ^ in the conduction band. In 
this way, 8 states are generated which can be expressed 
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in terms of the strong coupling eigenstates Eqs. ([6])- ([9]), 
see Appendix. It follows that in the strong coupling 
eigenbasis Eqs. the 2CK Hamiltonian |T]) takes 

the form of a generalized 2-channel Anderson impurity 
model, Eq. (|A.4[) . where the |\&°) play the role of the 
occupied and the |\P S ma) the role of the unoccupied im- 
purity. By a straight-forward Schrieffer- Wolff transor- 
mation^ for low energies, u> <C J, this Hamiltonian is 
projected onto the 2CK model 



H. 



(sc) 
2CK 



t 



E 



j 



C ■ C ■ A 

iacr jaa 1 2 



(10) 

where S = J2tt' F}&tt<F t , is the spin operator of the 
strong coupling compound, and J = (l/j)(4t) 2 /J, with 
7 = 30/46 1=3 2/3, is the effective spin flip coupling in the 
strong coupling regime (see the Appendix for a detailed 
derivation). Using, like in our NRG calculation of the 
following section, a flat DOS of po = 1 /4£, the dimen- 
sionless coupling reads, g = PqJ. Following Eq. the 
Kondo temperature is consequently given in the strong 
coupling regime by, 



(11) 



Comparison of Eq. (fTTjl with Eq. (|2|) exhibits the duality 
of the 2CK model in the weak and strong coupling limits 
via the identification 



1 

PoJ 



(12) 




FIG. 1: (Color online) Lowest energy levels of the isotropic 
2CK model as a function of the number of NRG iterations TV 
(N even) for different initial couplings J and A = 4 with 
900 states kept at each iteration. As a guide to the eye, 
the levels obtained in the different NRG iterations N are 
connected by straight lines. Independently of initial weak- 
coupling (J = 0.2), intermediate coupling (J = 0.7) or strong- 
coupling (J = 20) strengths, the same fixed point spectrum 
is reached. For odd numer of iterations N a non-equidistant 
fixed point spectrum is obtained as well (not shown). 
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Q 
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III. NRG TREATMENT AND RESULTS 

For the numerical solution of the 2CK problem we de- 
veloped an efficient NRG code, following Wilsons's orig- 
inal algorithm.— Since for the two-channel model the 
Hilbert space dimension grows particularly fast with the 
number N of NRG iterations, i.e. as 16 N , the use of 
conservation laws is essentual to reduce the Hamiltonian 
to block structure. The M-channel spin-i Kondo model 
has a full symmetry group of SU(2) spin x Sp(M), where 
Sp(M) is the symplectic group^ In the two-channel case 
(M = 2), the only decompositions into invariant sub- 
groups of Sp(2) are (i) SU(2) x U(l), corresponding to 
channel and charge conservation, and (ii) SU(2) x SU (2), 
corresponding to a separate axial charge conservation, 
used in the work of Pang and Cox4^ In our implemen- 
tation of the NRG for the 2CK model, we have chosen 
to use the decomposition (i), where we use the charge 
Q, z-component of the total (Kondo) spin S* ot and the 
z-component of the channel spin as labels for the many 
particle states only, corresponding to the following con- 



Thus we exploit only the U(l) subgroups of the full 
SU(2) spin and channel symmetries, respectively. Ac- 
cordingly, our code effectively uses a U(l) X U(l) X U(l) 
symmetry. This turned out to be an optimal compro- 
mize between computatinal efficiency and programming 
clarity. The Hamiltonians are diagonalized in each irre- 
ducible subspace \Q, S% ot , S* h ) and about 900 states were 
sufficient to be retained at each NRG iteration. After 
each NRG iteration the Hamiltonian is rescaled by the 
parameter y/A, A > 1^ The correct convergence of the 
NRG procedure was checked by comparing the results 
obtained with two different A-values, A = 3 and A = 4. 
It yielded excellent quantitative agreement, as seen below 
in Fig. [3 

We have solved the isotropic 2CK model for a wide 
range of bare spin couplings J in order to determine 
the J-dependence of Tr-. Typical flow diagrams of 
the energy eigenvalues are shown in Fig. [TJ exhibiting 
non-equidistant level spacings characteristic for the non- 
Fermi liquid fixed point j^A The fixed point coupling J* is 
characterized by the fact that, when the inital coupling is 
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FIG. 2: (Color online) Dependence of the Kondo temperature 
Tk on the bare coupling strength J (both in units of D), as 
determined by NRG for A = 3 and A = 4. The inset shows 
Tk on a logarithmic scale. 



J = J* , the energy eigenvalues settle immediately (after 
1 or 2 iterations) to their fixed point values (red dashed 
curves in Fig. [J), ft is thus identified as J* « 0.7-D in 
agreement with Ref. HH Following standard procedures, 
the Kondo temperature T% can be determined as the en- 
ergy scale either where the energy flow diagrams have 
an inflection point or where the first excited energy level 
has reached its fixed point value within, e.g., fO percent. 
Both definitions of this crossover scale are equivalent up 
to a constant prefactor, as seen for the weak coupling 
region in Fig. [3] Since, however, in the strong coupling 
region, J > J* , the complexity of the level flow makes it 
difficult to identify a single inflection point (see FigJTJ, 
we adopt the second definition. 

Our results for the dependence of Tk on the bare 
Kondo coupling J are shown in Fig. [5] ft shows a strong 
peak at around J = 0.7-D, as expected. The results for 
the two discretizations considered, A = 3, A = 4, show no 
significant differences. The deviations in the intermedi- 
ate coupling regime around the peak maximum in Fig. [2] 
arise from the difficulty to determine the exact T K when 
the crossover happens at the very beginning of the NRG 
iterations where the energy resolution is low. The behav- 
ior of T K is examined over nearly three decades of J and 
extends over more than fO decades in T K , as illustrated 
in the inset of Fig. O 

The J-dependence of Tk can be further analysed by 
plotting it in Fig. [3] versus the parameter — (l/<7 + 7s). ft 
shows the exponential behavior of Tk as a function of 1/ J 
in the weak coupling limit 1/ J — > oo and as a function of 
J in the strong coupling limit I /J — > oo, with logarithmic 
corrections towards the intermediate coupling regime, in 
agreement with Eqs.[2land[TTl Note that the strong cou- 
pling and the weak coupling branches in Fig. [3] are par- 
allel to each other, i.e. the NRG quantitaively confirms 
the analytical value of the parameter 7 = 30/46. The 




_l 1 1 1 1 I 1 1 1 1 I 1 1 1 1 L_ 

-20 -15 -10 -5 

-(1/9+79) 



FIG. 3: (Color online) The Kondo temperature Tk is shown 
on a logarithmic scale versus the parameter —(^7 + JpoJ), 
7 = 30/46 sa 2/3. The upper branch of the curves corresponds 
to the weak coupling, the lower branch to the strong coupling 
regime. The results for Tk obtained from the "inflection point 
method" (see text) in the weak coupling regime (J < J*) are 
shown for comparison and differ only by a constant prefactor. 



shift of the two branches can be traced back to the fact 
the strong coupling Anderson impurity model, Eq. (|A.4p . 
produces higher-order logarithmic corrections which are 
different from those of the weak coupling model, Eq. (JTJ) 
and which are, thus, not included in the effective low- 
energy 2CK model, Eq. (JTDJ). 



IV. CONCLUSION 

The two-channel Kondo model exhibits for low ener- 
gies a duality between the regions of weak and strong 
bare Kondo coupling J. This results from the fact that 
in both limits, J — > and J — > 00 the ground state is 
doubly degenerate. While for J — > it is the decou- 
pled impurity spin doublet, for J — * 00 it is a doubly 
degenerate quantum frustrated 3-body state, comprised 
of the impurity spin and the conduction electron spins 
located at the impurity site in each of the two channels. 
We have shown that, hence, the complete strong coupling 
behavior can be obtained from the solution in the weak 
coupling regime via the identification of the dimension- 
less coupling, jg — > 1/g, where 7 = 30/46 w 2/3. These 
results have been confirmed quantitatively by the exact 
numerical renormalization group solution of the problem. 

As a result, the dependence of the Kondo tempera- 
ture Tk on the J is strongly peaked at the two-channel 
Kondo fixed point coupling, J = J* « 0.7, and decays 
exponentially both for small and for large couplings. The 
maximum is of the order of the band cutoff, Tk(J* ) ~ D, 
with non- Fermi liquid behavior for all energies below Tr. 

We conjecture that this could be the reason why in ex- 
perimental conductance anomalies of nanoconstrictions 
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with two-channel Kondo signature a 17 ' 18 i 21 no broad dis- 
tribution of Tx is observed: The band cutoff and hence 
Tk{J*) hi two-channel Kondo systems can be provided 
by a decoherence scale of the order of a few Kelvin^ This 
would mean that, even if there is a broad distribution of 
bare couplings, only for those couplings sufficiently close 
to J* the non-Fermi liquid behavior would extend up 
to sufficiently high energies to be observable. However, 
more detailed calculations as well as a detailed micro- 
scopic model for the two-channel Kondo defects will be 
required to substantiate this conjecture. 
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APPENDIX: DETAILS ON THE DUALITY 
ANALYSIS 
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Destruction of an electron from the 3-particle com- 
pound ground states (Eqs. Q}, (O in channel a = ±1 
generates the (unnormalized) states, 



co T i|*?> = 
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which can be expressed in terms of the strong coupling singlet/triplet eigenstates Eqs. 
bosonic creation operators for the latter states, 



I*. 



B] m& \vac) 



as indicated. We define 
(A.2) 



which transform with respect to the channel SU(2) group according to the adjoint representation, i.e. a = —a. 
Together with the fermionic operators of Eqs. Q, (O they satisfy the constraint 



/ , a a ' / j sma J -' sma x i 



(A.3) 



an expression of the uniqueness of the strong coupling basis states. In the strong coupling basis, using Eq. (|A.1[) the 
2CK Hamiltonian {!]) then takes form of a generalized two-channel Anderson impurity model in one dimension, 



H = 1 E £ c L<w- j Y, f ° f °--a j Y. b 1^ b ^ + \ j E B l 



ma Ima 
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z=±l a 
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2 Cial^ooa^i ^g t mt- D i-ia- r | 
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-t Rt 



. ial B{ m F l+ H.c. 



(A.4) 



where V = 2t plays the role of the band-impurity hybridization and the factor 2 arises from the fact that there 
is hopping from the impurity site to the two sites i = ±1. By means of a Schrieffer- Wolff transformation 32 this 
Hamiltonian maps for energies u> <C J onto the effective 2CK model Eq. (TTU|) . where potential scattering terms have 
been neglected. Since only the intermediate (bosonic) states with m = contribute to an effective Kondo spin flip 
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(i.e. only products of the 1st and the 4th term and of the 
the effective spin flip coupling, as defined through Eq. (TT 




where 7 = 30/46 « 2/3. 



d and the 6th term of the hybridization part in Eq. (IA.4I) ~1 . 
reads, 
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